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Abstract: We solve the Dyson-Schwinger equations of the ghost and gluon propagators of Landau 
gauge Yang-Mills theory together with that of the ghost-gluon vertex. The latter plays a central role 
in many truncation schemes for functional equations. By including it dynamically we can determine 
its influence on the propagators. We also suggest a new model for the three-gluon vertex motivated 
by lattice data which plays a crucial role to obtain stable solutions when the ghost-gluon vertex is 
included. We find that both vertices have a sizable quantitative impact on the mid-momentum regime 
and contribute to the reduction of the gap between lattice and Dyson-Schwinger equation results. 
Furthermore, we establish that the three-gluon vertex dressing turns negative at low momenta as 
suggested by lattice results in three dimensions. 
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1 Introduction 



The propagators of Yang-Mills theory have been widely investigated because, on one hand, they 
provide the starting point for calculations of hadronic observables and investigations of the phase 
diagram of strongly interacting matter as described by quantum chromodynamics (QCD). On the 
other hand, they provide some direct insight into the mechanisms behind non-perturbative phenomena 
like confinement and dynamical chiral symmetry breaking. Typically, in the Landau gauge the ghost- 
gluon vertex has been a pivotal object in such studies when using functional methods, since its simple 
structure is the basis for most truncation schemes. The original motivation for this was Taylor's 
non-renormalization argument for the Landau gauge [1]. The employed truncation schemes were so 
successful that the complications arising by including the vertex dynamically into numeric calculations 
have been postponed for some time. In the meantime additional information was gathered [2-9] with 
functional methods backed up by calculations on the lattice [10-13] that solidified the reliability of 
the used truncation schemes and supported the expectations that the modifications induced by a self- 
consistent inclusion of the vertex are small and only on the quantitative level. Here we test these 
expectations by including the ghost-gluon vertex dynamically into the system of equations we solve, 
in order to assess any quantitative or even qualitative differences. 

While major qualitative changes are certainly not expected, the extent to which quantitative 
modifications occur is not clear. Since the ghost-gluon vertex is an important part of the Yang-Mills 
sector, we should study it as well as possible to ensure that we fully understand the non-perturbative 
infrared regime of Landau gauge QCD. Furthermore, it recently turned out that a simple ghost-gluon 
vertex as employed in the vacuum is insufficient for non-zero temperature calculations [14]. Finally, 
it is desirable to improve existing results for the propagators quantitatively in order to become more 
competitive with lattice results. The inclusion of the ghost-gluon vertex constitutes one step in this 
direction, also because it is a prerequisite for the investigation of other quantities like the gluonic 
vertices. The lowest one, the three-gluon vertex, is expected to have even more quantitative influence 
on the propagators. In order to test this assumption, we investigate a new three-gluon vertex model. It 
is Bose symmetric, has the correct ultraviolet (UV) anomalous dimension and features a zero crossing 
as observed in three-dimensional lattice data. The latter property has not been clearly observed in 
four dimensions yet. Thus we also compute the infrared (IR) leading contribution of the three-gluon 
vertex for one momentum configuration to confirm this. 

Perturbatively the ghost-gluon vertex is well studied. A detailed one-loop calculation for general 
external momenta can be found in Ref. [15] and a three-loop calculation for the so-called asymmetric 
point in Ref. [16]. An often used argument about its IR behavior goes back to Taylor [1]. We will 
comment on it in more detail below in Sec. 2.2. Further analytic studies have been performed with 
Dyson-Schwinger equations (DSEs) [5, 6] and first functional renormalization group calculations for 
non-zero temperature were also done [14]. On the lattice the vertex was investigated in Refs. [10-13], 
but unfortunately the results are by far not as detailed as for the propagators. Still, they are the 
most reliable source of quantitative information. Semi-perturbative results from DSEs can be found in 
[3, 17]. Here we will present the first calculation with full momentum dependence and back coupling 
effects. 

The dynamical inclusion of the ghost-gluon vertex required at the same time to improve on the 
three-gluon model vertex. The reason is that in the gluon propagator DSE the contributions of the 
ghost and gluon loops naturally balance to keep the gluon dressing function positive. If we now 
change the ghost-gluon vertex this balance might get lost and the gluon loop could take over in the 
mid-momentum regime rendering the gluon propagator negative. This gives a natural constraint on 
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any vertex model and was one motivation for introducing the new model for the three-gluon vertex 
that takes into account the known properties from lattice [13] and Dyson-Schwinger studies [6, 18]. 
The improved vertex has a considerable effect in the mid-momentum regime, where the difference to 
lattice calculations is most evident. 

We start with an overview of the Green functions employed here and their DSEs in Sec. 2. There 
also some remarks about the Taylor theorem can be found and the model for the three-gluon vertex 
is explained. The employed methods are reviewed in Sec. 3. Results are presented in Sees. 4 and 5. 
We conclude with a summary in Sec. 6. In two appendices we present the kernels of the DSEs and 
further details on the UV behavior of the propagators. 

2 The system of Dyson-Schwinger equations 
2.1 The two-point Dyson-Schwinger equations 

The DSEs for the ghost and gluon dressing functions, G(p 2 ) and Z(p 2 ), respectively, read 

G ^y = Z 3 + N c g 2 Z 1 J^Z(q 2 )G((p + q) 2 )K G (p, q)r A£c (q-p + q,p) (2.1) 



' Z- i + N c g 2 ~Z 1 I G(q 2 )G({p + q) 2 )K 9 z h (p, q)T Asc (p;p + q, q) 



Z{p 2 ) 

+ N c g 2 Z 1 I Z{q 2 )Z{(p + q) 2 )K%{p,q)T A \p,q,p + q), (2.2) 



'/ 



where J stands for J d q/(2ir) and Z\ and Z\ are the renormalization constants of the ghost-gluon 

and three-gluon vertices, respectively, whereas Z 3 is the renormalization constant of the ghost propa- 
gator and Z 3 that for the gluon propagator. The quantities T Acc (k;p, q) and T A (k,p,q) are dressing 
functions of the ghost-gluon and three-gluon vertices, respectively. The two DSEs are shown dia- 
grammatically in Fig. 1. The ghost and gluon propagators D g i l (p 2 ) and D g i tllv {p 2 ) and their dressing 
functions are related by 

D^p 2 )-- ^!, D gl ^(p 2 ) := P (p)?^l. (2.3) 
p l p l 

The kernels Kq, K^ h and K% are given explicitly in Appendix A. The coupled system for the two 
propagators was investigated many times in the literature with functional methods, see, for example, 
[2, 5, 7, 19-31]. Recently it has even been solved in the complex p 2 -plane to extract the corresponding 
spectral functions for gluons and ghosts in complete agreement with local quantum field theory but 
with the expected positivity violations [32]. For complementary Euclidean results from Monte Carlo 
simulations see, e.g., [33-46]. 

In Section 5 we will compare our results also to lattice data. Thus we shortly want to recall a 
few facts concerning the calculation of propagators on the lattice. It is often claimed that all lattice 
results show an unequivocal picture for the propagators. However, several open issues exist that need 
clarification. For example, it is well known that Gribov copies have an effect on the propagators 
[12, 34, 47-49]. An interpretation of this in terms of non-perturbative gauge completion was put 
forward in Ref. [48]. With functional equations one finds two types of solutions [5, 28, 50] characterized 
by the value of the ghost dressing function at zero momentum: If it is finite, the family of decoupling 
solutions emerges which has a non-zero gluon propagator at zero momentum [5, 27, 28, 50, 51]. The 
second type is the scaling solution, for which the ghost dressing diverges and the gluon propagator 
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Figure 1. Full two-point Dyson-Schwinger equations of Landau gauge Yang-Mills theory. All internal propa- 
gators are dressed. Thick blobs denote dressed vertices. Wiggly lines are gluons, dashed ones ghosts. 



vanishes [19]. Of course it is tempting to look for a direct correspondence between the choice of Gribov 
copies and the family of solutions. However, Monte Carlo methods sample different Gribov copies in a 
way not accessible to functional methods, so that unveiling such a relation is in no way trivial. Recently 
promising new results in this direction were presented [52]: Choosing Gribov copies by the value of 
the lowest eigenvalue of the Faddeev-Popov operator seems to lead to changes in the propagators 
akin to those obtained when varying the boundary conditions in the DSEs. Another point in need 
of better understanding is the behavior in the strong coupling regime, for which all momenta can be 
considered below Aqcd and thus in the IR. There two solution branches are found with characteristics 
of decoupling and scaling solutions [53-55]. Due to these ambiguities it is in principle important to 
consider the details of the gauge-fixing algorithm used to sample the Gribov copies of Landau gauge 
on the lattice when comparing the infrared behavior of the propagators with lattice data. All lattice 
data used for comparisons herein were obtained from sampling, with some algorithm-specific bias, 
local minima of the standard gauge-fixing potential for lattice Landau gauge thus yielding decoupling. 

2.2 The ghost-gluon vertex Dyson-Schwinger equations 

The ghost-gluon vertex has two different DSEs which differ by the field that is attached to the bare 
vertex, see Fig. 2. Of course both full equations are equivalent and we can choose the one most suitable 
for our calculations. It seems tempting to employ the DSE with the ghost leg attached to the bare 
vertex, because it has only four diagrams. However, both UV and IR consistent truncations depend 
on the dressed three-gluon vertex for which we use a model. To make our calculations as much model 
independent as possible we thus want to avoid the dressed three-gluon vertex in the ghost-gluon vertex 
DSE. We therefore consider the first DSE in Fig. 2. 

In the UV the leading diagrams are the two triangle diagrams. The two other one-loop diagrams, 
the gluon and the ghost loops with quartic ghost-gluon and ghost interactions, respectively, are UV 
subleading because the four-point functions therein are not primitively divergent. At the same time 
the triangle diagrams contain the leading IR corrections. This can be seen most easily for the scaling 
type solution: There, the two IR leading diagrams as indicated by the power counting scheme [18] are 
the two diagrams with a bare ghost-gluon vertex. However, one of them contains a quartic ghost vertex 
which behaves as (p 2 ) K in the IR. We expect the contribution from this diagram to be suppressed as 
compared to the triangle diagram with a bare ghost-gluon vertex because there is no tree-level four- 
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Figure 2. The ghost-gluon vertex DSEs. All internal propagators are dressed. Top: Version with the bare 
vertices attached to the gluon leg. Bottom: Version with bare vertices attached to the ghost leg. Thick blobs 
denote dressed vertices. Wiggly lines are gluons, dotted ones ghosts. 



ghost coupling in Landau gauge. It corresponds to a two-loop contribution as can be seen explicitly by 
inserting the DSE of the quartic ghost vertex. For the decoupling type solution no IR enhancements of 
the vertices are expected [5], so the main corrections should also come from the one-loop terms. Note 
that this truncation is almost the same as the first three diagrams of the other ghost-gluon vertex 
DSE depicted at the bottom of Fig. 2 except for the dressed three-gluon vertex. Consequently all 
information of the diagrams neglected in our truncation is contained in the full three-gluon vertex and 
the ghost-gluon scattering kernel. 

The dressed ghost-gluon vertex is described by two dressing functions. A typical choice is 

rf c ' abc (fc; P ,g) := f abc Tf c (k;p,q) := 19 f abc (p„A(k;p,q) + k^B{k-p,q)) . (2.4) 

Here the momenta k, p, q refer to the gluon, anti-ghost and ghost legs, respectively, as indicated in the 
superscript of Y Acc . By convention we have all momenta ingoing. The dressing A(k;p, q) is symmetric 
under exchange of the ghost momenta, A(k,p, q) = A(k;q,p), but B{k;p,q) is not [2]. In the Landau 
gauge only A(k;p, q) contributes to Green functions due to the transversality of the gluon propagator. 
However, it is slightly misleading to call A(k;p,q) the transverse dressing of the ghost-gluon vertex, 
as is often done. A clear separation between transverse and longitudinal parts can be achieved with 
[2, 56]: 

rf c > abc (k;p, q) :=tgf abc (P^(k)p„Df c (k;p, q) + k^Df c %k;p, q)) , (2.5) 
where P M1 ,(fc) is the transverse projector g^u — k^k^/k 2 . The relations between the two sets of dressings 
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are 



Df c (k;p,q)=A(k;p,q), (2.6) 
Df c (k;p, q) = B(k;p, q) + A(k;p, q)p- k/k 2 (2.7) 

from which one sees that A(k;p, q) also appears in the longitudinal dressing function and the transverse 
part is symmetric under exchange of the ghost momenta, Df cc (k;p, q) = Df cc (k; q,p). The bare vertex 
is obtained for A(k;p, q) = 1, B(k;p, q) = and Df cc (k- 1 p 1 q) — 1, Df cc (k;p, q) = k-p/k 2 , respectively. 
Note that the choice of a ghost-gluon vertex that ensures the IR transversality of the gluon DSE as in 
Refs. [2, 28] corresponds to setting the longitudinal dressing Df cc (k;p,q) to zero in the IR. However, 
when using the transversely projected gluon DSE Df cc (k;p, q) does not enter anyway and only appears 
when testing if the longitudinal components of the gluon propagator vanish. 

Using the transverse/longitudinal basis one can easily see that the ghost-gluon vertex has no 
transverse part for zero ghost momentum: 

rf c (-p;p,0) =ig(P^(-p)p u Df cc (-p-,p,0)-p fl Df- cc (-p-,p,0)) =-igp M Df Sc (-p;p,0). (2.8) 

Or put in terms of A and B: 

rf c (-p;p,0) = ig f ahc p, (A(-p;p,Q) - B(-p;p,0)) . (2.9) 

This observation is related to what has become known as Taylor's non-renormalization argument [1]. 
Taylor argued that due to transversality one can write the ghost-gluon vertex as 

Tf c (k,p,q) = igp^+p v q p F^ p {k,p,q). (2.10) 

For p 2 = q 2 = k 2 = the second term should equal zero because of the factor q p [57] leading to 
Z\ — 1 if the vertex is rcnormalized by momentum subtraction at p 2 = q 2 = k 2 = 0. 1 However, this 
statement implicitly assumes that the value of the ghost-gluon vertex is the same independent from 
which direction the IR is approached. For the scaling solution it was already observed in Refs. [3, 4, 6] 
that the limits p 2 , q 2 , k 2 — > and q 2 — > 0, p 2 — > are not equal. Furthermore, any possible singularities 
in F^ p ^,q,p) ar e disregarded. However, since there are no indications of this, we ignore this issue in 
the following. Then the ghost-gluon vertex becomes bare in the limit of vanishing incoming ghost 
momentum q. More directly this can be inferred from the ghost-gluon vertex Slavnov-Taylor identity 
(STI) [20]: 

-^y(Al(y)c c (z)c a (x)) - -^(Al{z)c\y)c a (x)) = \g f ade (c d (xy(x)c?(z)<?(y)). (2.11) 

If one ignores the connected part of the ghost scattering kernel one arrives at 

ik li Y^k-p,q)G- 1 {k 2 )+iq li Y^q-p,k)G- 1 {q 2 ) = gp 2 G- 1 (p 2 ). (2.12) 

If the connected part is maintained, the following conclusions remain consistent but a rigorous proof 
is lacking. Of course, since this is an STI, it only provides information about the longitudinal part of 
the ghost-gluon vertex. Assuming that the longitudinal dressing function Df cc (k;p,q) has no strong 

1 Note that Z\ = 1 is automatically true for renormalization schemes based on minimal subtraction like MS [58] or 
MiniMOM [59]. 
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divergence for k — > 0, we can take q or k to zero and use lim 9 2_>,0 Q 2 G 1 (q 2 ) = 0, which is valid for 
both decoupling and scaling solutions, to obtain 

Df c (-p;p,0) = -l. (2.13) 

This is equivalent to 

r„(-p;p,0)=igp li , (2.14) 

i.e., the ghost-gluon vertex reduces to the tree-level contribution in this limit as expected. From 
Eq. (2.13) we can switch to the dressings A(k;p, q) and B(k;p,q) and obtain 

B(-p;p,0) - A(-p;p,0) = -1. (2.15) 

This relation was noted also, for example, in Ref. [58]. Interestingly, it also means that the transverse 
dressing is not necessarily zero, because of Df cc (— p;p, 0) = A(— p;p, 0). Thus when we calculate 
A(k;p, q) or Df cc (k;p, q) it will not vanish for q — > and the vanishing of the radiative corrections to 
the ghost-gluon vertex for zero ghost momentum is only of kinematic origin. 

In order to calculate the ghost-gluon vertex one can either project out the dressings to obtain 
scalar equations, which is most easily done for the Dffi c basis, or decompose the appearing vector 
integrals into scalar ones, which is the method of choice for the A/B basis. Here we follow the latter 
approach. The truncated DSE for the vertex reads then 

A{k\P,q) = 
1 , N <9 2 



A(p, k) J r 
N c g 2 
A(p,k) 



K? cc (k,p, r)D gh (r 2 )D gh ((r + k) 2 )D gl {(p - r) 2 )A{r - p; p, -r)A(p - r; k + r, -k - p) 



K^ c (k,p,r)D gh ((r-p) 2 )D gl (r 2 )D gl ((r + k) 2 )A(r;p,p- r)A{k + r;p-r,-k- p). 

(2.16) 



The kernels K^ cc (k,p,r) and K2 CC (k,p 7 r) are given in Appendix A. Since the dressing B(k;p,q) is 
not required for the Landau gauge, we do not give the corresponding expressions. 

2.3 Three-gluon vertex 

For the three-gluon vertex we will employ an ansatz. In the past different expressions have been used. 
The main guideline in their construction was the correct UV behavior for the gluon propagator, but 
also other aspects like quadratic UV divergences in the gluon DSE played a role, see Sec. 3.1 and, for 
instance, Refs. [20, 24, 28, 60, 61]. Since we use the transversely projected gluon two-point DSE we 
only deal with the completely transverse part of the vertex, which is unrestricted by Slavnov- Taylor 
identities. Quantitatively the vertex is expected to contribute to the closing of the gap between DSE 
and lattice results for the gluon dressing function in the mid-momentum regime. Another source of 
missing strength there are two-loop diagrams, of which one, the so-called squint diagram (the last 
one in Fig. 1) also contains the three-gluon vertex. Using the functional renormalization group the 
propagators can be calculated such that this gap practically vanishes [28], but due to the different 
structure of the corresponding equations, it can not be inferred how much additional support in the 
mid-momentum regime is expected from which source in DSEs: two-loop diagrams or the dressing 
of the three-gluon vertex in the gluon loop. First calculations of two-loop diagrams in DSEs were 
performed in [62, 63]. 
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The IR behavior of the three- gluon vertex shows some interesting features. For the scaling type 
solution it was found that the uniform IR limit, i.e., when all momenta approach zero simultaneously, 
obeys the power law (p 2 )~ 3K [18] and that possible kinematic soft divergences can occur that go like 
1 — 2k [5]. However, for the transversely projected vertex the degree of divergence is reduced to 
3/2 — 2k [6, 7] and such divergences play no role. For decoupling type solutions no IR divergences in 
the three-gluon vertex are expected [5]. 

In light of the discussion of possible ambiguities in the IR behavior of the solutions of Landau 
gauge Yang-Mills theory it is interesting to consider the three-gluon vertex as obtained with Monte 
Carlo simulations [13]: While in four dimensions the statistics is not high enough to allow a detailed 
assessment of the IR behavior, the data in three dimensions shows a zero crossing and it becomes 
negative. Also in two dimensions [64] a zero crossing is found and additionally the data agrees with 
the expected power law from a scaling solution [56, 65]. This is not unexpected, because it has become 
by now clear that in two dimensions no decoupling solution exists [56, 66, 67]. However, in three and 
four dimensions only the decoupling type solution is found on the lattice and hence the realization of 
such a power law would indeed be unexpected. Since no IR singularities are expected [5] , the vertex 
most likely approaches a constant but negative value at zero momenta. Further lattice results are 
therefore desirable for clarification of the exact IR behavior of the vertex. Based on the qualitative 
similarity of the propagators of three and four dimensions we assume that the zero crossing also exists 
in four dimensions and it is therefore included in our vertex model here as well. Further evidence in 
support of this is presented in Sec. 4 from a DSE calculation of the three-gluon vertex. 

The employed construction of the three-gluon vertex is based on the one introduced in Refs. 
[24, 60, 61]: 

D A (p, q,-p-q) = — rr , (2.17) 

P q> Z x [Z(y)Z(z)] 1+a ' V ' 

where Z\ is the renormalization constant of the three-gluon vertex, y — q 2 , z = (p+q) 2 and 5 — —9/44 
is the anomalous ghost dimension, related to that of the gluon by 1 + 7 + 25 = 0. a is a parameter 
which is irrelevant in the UV where it drops out. A typical choice is a = 35, which renders the 
dressing finite at vanishing momenta in case of the scaling solution [60]. For the decoupling type 
solution the corresponding choice is a — — 1. For our construction described below two observations 
are important: First, this vertex construction is not Bose symmetric; the momenta p and q are the 
external and internal momenta, respectively, in the gluon loop. Secondly, it contains the inverse 
of the renormalization constant of the three-gluon vertex to cancel the Z\ appearing in the gluon 
loop of the gluon DSE due to the bare three-gluon vertex. By construction this model ensures the 
corrected logarithmic running of the gluon loop according to resummed perturbation theory. However, 
it also means that the anomalous dimension of the three-gluon vertex model itself is off by a factor 
of two. In earlier works this was remedied by replacing Z\ by a momentum dependent function 
[20]. Because the correct anomalous dimension requires also contributions of order g 4 and higher, 
this cannot be reproduced directly with the employed truncation of the gluon DSE and we will use a 
similar renormalization group improvement term. Other models with Bose symmetry and the correct 
anomalous dimension of the vertex were used in Ref. [68], in the context of the quark-gluon vertex 
DSE, and in Ref. [30] for the Yang-Mills propagator system. However, there the anomalous dimension 
of the gluon propagator was only reproduced by interpolating between these models in the IR and the 
one of Eq. (2.17) with a — in the UV [69] so that in effect the employed model did not have the 
correct anomalous dimension. 



- 8 - 



We first address the issue of Bose symmetry. Naively adding factors of G(x) and Z(x) with x = p 2 
does not work, but instead we use a Bose symmetric combination of momenta: 2 

D^(p, q,-p - q) = G (^f^) ° Z (^f^) " ■ (2.18) 

In order to determine a and /3 we demand that the vertex has the correct anomalous dimension: 

5a + jf3 = l3g . (2.19) 

From the Slavnov- Taylor identity Z\ = Z-^jZ^ we know 73 9 = 1 + 35 = 17/44. Furthermore we require 
IR finiteness, i.e., 

5 C a + S A P = 0, (2.20) 

where 5 C and 5a are the IR exponents of the ghost and gluon propagators, respectively. This leads 
for the scaling solution to a = —2 — 65 and f3 = —1—35. For decoupling we have a = 3 + 1/5 and 
j3 = 0. Note that IR finiteness is not per se required, since we will control the IR behavior by another 
term, but for the current model it is advantageous that D A ' UV (p, q, —p — q) does not vanish at zero 
momentum. 

Next, we modify the IR behavior in order to incorporate the zero crossing expected from lattice 
calculations by the following expression: 

D A3 ' IR (p, q, -p-q) = h IR G(x + y + zf{f 9 (x)f° {y)f 9 (z))\ (2.21) 

with the damping factors 

A 2 

fa{x) : = 39 . (2.22) 

A similar construction to obtain the zero crossing was employed in two dimensions [56] . The exponent 
of the damping functions is 4 in order to reproduce the fast change seen in three-dimensional lattice 
calculations. The impact of this IR part is controlled by the two parameters hm and A3 g . Naturally 
hm is negative to produce a zero crossing. From currently available lattice data [13] one can estimate 
bounds for the parameters, which, however, is not a straightforward procedure because the conversion 
from internal to physical units can only be done a posteriori. The total vertex employed in this study 
is given by 

D A \ Pl q, -p-q) = D A3 > IR (p, q, -p - q) + D A ^ UV (p, q, -p - q). (2.23) 

Finally, we want to incorporate the correct UV behavior of the gluon propagator by adding a 
corresponding improvement term in the gluon loop. As mentioned above, this is necessary because 
the bare vertex in the gluon loop does not give any contribution to the anomalous dimension. Note 
that this problem would be absent, if we had both three-gluon vertices dressed. This is the case for 
the functional rcnormalization group [28] and also for equations of motion obtained from nPI effective 
actions [70] . Earlier studies accounted for this in the construction of the three-gluon vertex model by 
using an anomalous dimension a factor two too large [24, 28, 60, 61]. However, since we want to use 
lattice calculations as input to improve on the vertex model, we have to split the renormalization group 



We acknowledge discussions with Christian S. Fischer and Stefan Strauss who employ a similar construction. 
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improvement term and the three-gluon vertex model. In analogy to the three-gluon vertex model we 
choose a Bose symmetric expression: 

d£g(p, 1, -P-Q) = ^D A3 > uv (p, q, -p-q), (2.24) 

This term and the three-gluon vertex model are combined in the function T a3 (p, q, —p — q) in the 
gluon DSE: 

T A " (p, q,-p-q)=D A \p,q,-p-q)Df G (p, qi -p- q ) (2.25) 

Here it becomes obvious why we required IR finiteness of the UV part, because otherwise the whole 
expression could vanish for zero momenta. The presented model for the three-gluon vertex is compared 
with lattice data in Section 4. 

3 Methods 

3.1 Quadratic divergences in the gluon DSE 

The use of a UV regularization which is not gauge invariant entails that the gluon DSE is plagued by 
quadratic divergences. Several ways of dealing with them have been employed in the literature. For 
example, subtractions at the level of the integrands [60, 61] or modifications of the vertices [28] were 
used. Another possibility is the projection with the so-called Brown-Pennington projector [30, 69, 71], 
which, however, leads to the dependence on the longitudinal part of the ghost-gluon vertex. It also 
introduces a spurious dependence on an additional gauge parameter [2, 72], which should be irrelevant 
in the Landau gauge. Thus we employ the usual transverse projection. 

Here we use the following method to handle the quadratic divergences: Instead of adding an 
additional term in the kernel of the gluon loop to subtract the quadratic divergences of both loops as 
done, for instance, in Refs. [60, 61], we add such terms for each loop separately. This has, for example, 
the advantage that truncations that take only the ghost loop into account can be used. Another benefit 
is that a standard fixed-point iteration becomes a viable solution method. If the quadratic divergences 
are subtracted only via the gluon loop, instabilities are introduced in the iterative procedure and 
no solution is found [73]. This made the Newton method a crucial aspect of most previous solution 
processes. 

Since quadratic divergences are a problem of the UV regime, we want to avoid any interference 
with the IR part. Thus the additional terms are suppressed in the IR by a damping factor. The 
corresponding expressions can be found in Appendix B. 

3.2 Solution method 

The system of the three DSEs is split into two subsystems, which are iterated in turn until a given 
precision is reached. The first subsystem consists of the propagator DSEs, for which we employ the 
Newton method. This method is used throughout the literature and we refer for details, for example, 
to Refs. [21, 24, 60, 73, 74]. Alternatively we used a direct fixed point iteration for all three equations. 
We solved the equations on CPUs, but note that recently it was demonstrated that implementing such 
a system on GPUs can lead to a significant speedup [75]. 

For the ghost-gluon vertex a standard fixed point iteration is employed. Therefor we put the 
ghost-gluon dressing function A[k; q,p) on a grid and use linear interpolation for intermediate points. 
If points outside the grid are encountered we take the value at the boundary, since the vertex is 
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constant in the IR and the UV. The validity of this procedure can also be confirmed a posteriori from 
the results. As variables we choose the momenta squared of the gluon and the anti-ghost, k 2 and q 2 , 
respectively, and the cosine of the angle, c t = cos#, between these two momenta. The allowed intervals 
for the variables are then [0, oo], [0, oo] and [—1, 1]. If we used a basis with three momenta squared, 
the intervals would depend on the values of the other variables through momentum conservation which 
would complicate the calculations unnecessarily; hence this choice is here more advantageous. Since 
the dressing function of the ghost-gluon vertex has to be calculated on a three-dimensional grid, the 
computing time increases considerably when compared with the propagator system. Furthermore, the 
kernels are more complicated and the integrations are three- instead of two-dimensional, with two angle 
variables 6\ and 62- For the propagators in total three integrations have to be performed, because 
the radial momentum integration is split at the value of the external momentum. The integration of 
the vertex requires similar splittings [17] at r 2 — p 2 with #2 = and at r 2 — k 2 with Q\ — ir and 
62 = 7r — 6, where r is the radial integration variable and 9\ and 02 are the angle integration variables. 
Thus in total six integrations have to be done. 

The two subsystems are iterated in turn until the required precision is reached. The framework 
CrasyDSE [73] was employed for all calculations. Color factors were used throughout for the gauge 
group SU(S). 

4 Effects of the three-gluon vertex 

In this section we illustrate the quantitative effects of the three-gluon vertex model and investigate the 
impact of the position of its zero crossing by considering the system of propagators. The ghost-gluon 
vertex is set to tree-level. By construction the UV behavior of the propagators is not subject to any 
change. The most crucial modifications will be seen in the mid-momentum regime. For the IR only 
quantitative changes are expected, like a varying value of the gluon propagator at zero momentum. 

We start by comparing three different three-gluon vertex models: The one of Eq. (2.17), the Bose- 
symmetrized one (2.18), and the IR enhanced one, viz. Eq. (2.23) with a non-zero hm. The results for 
a decoupling type solution are shown in Fig. 3. As can clearly be seen, both the Bose symmetrization 
and the IR enhancement enlarge the bump in the gluon dressing function. Also the zero value of 
the gluon propagator changes and raises to higher values. The ghost propagator, on the other hand, 
remains nearly unaffected. In Section 5 we will compare the propagators also to lattice results and 
will see that the enlarged bump considerably reduces the gap between DSE and lattice results. 

What we need for the three-gluon vertex model are values for its parameters. Ideally we could 
get them from lattice data. However, due the presently large uncertainties this is still somewhat 
ambiguous. Thus we test several values and their influence on the propagators. For simplicity we fix 
the parameter hm = —1 and only vary A3 g . In Fig. 4 we show results for values of As g between and 
1.7 GeV. The choice A3 S = OGeV means that the IR enhanced part is dropped. In addition to the 
propagator dressing functions we also show a comparison of the three-gluon vertex model with lattice 
data at the symmetric point. This data is for SU(2), but since it is known that the propagators look 
extremely similar for SU(2) and SU(3) [76] one might expect the same for the three-gluon vertex. 
The two other available configurations (orthogonal momenta of equal magnitude, or one momentum 
vanishing) look similar. One can see from the lattice data that the turnover to negative values has to 
proceed rather rapidly. This is the reason why the exponent of the damping functions in Eq. (2.21) 
is chosen as 4. For the plot of the three-gluon vertex in Fig. 4 the final propagators were used in the 
model. As can be seen, the shown models have a zero crossing at higher values than indicated by the 
data. Such high values were chosen to demonstrate the effect it has on the gluon dressing function: 
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Figure 3. Comparison of different three-gluon vertex models: model of Eq. (2.17) (blue/dotted line), Bose 
symmetrized model (Eq. (2.23) with hiR — 0; green/dashed line), Bose symmetrized model with IR enhance- 
ment (Eq. (2.23) with hm — — 1, A3 9 = 1.7 GeV; red/continuous line). Top: Ghion dressing function (left) 
and propagator (right). Bottom: Ghost dressing function, for which the three lines almost lie on top of each 
other. 

Only for a zero crossing at momenta high enough the gluon dressing is affected. This is evident in 
the plot of the gluon dressing functions where the first three curves, corresponding to no zero crossing 
and crossings below 600 MeV, show not much difference. For higher values of the zero crossing the 
bump increases faster. Again we see no effect on the ghost dressing function and only a small one on 
the zero momentum value of the gluon propagator. 

The observations on the dependence on the place of the zero crossing lead to the conclusion that 
the actual value of the three-gluon vertex in the deep IR is not very important for the propagators. 
Indeed, for values of up to roughly 0.7 GeV for A 3ff we do not find any difference to the propagators 
obtained with a three-gluon vertex without zero crossing. This is not too surprising considering that 
the three-gluon vertex in the gluon loop is multiplied by two gluon propagators. For the scaling 
solution this leads to an IR suppression of (p 2 ) K and for decoupling solutions to (p 2 ) 2 ■ 

In order to shed more light on the question of the existence of a zero crossing, we calculated the 
IR leading part of the three-gluon vertex for the symmetric momentum configuration. In the IR the 
so-called ghost triangle is the dominant contribution. Since it contains no three-gluon vertices there 
are no back coupling effects and the calculation consists of only one iteration. Of course, the mid- 
momentum and UV behavior are not described accurately in this rough truncation. Using the ghost 
propagators from the propagator calculations and bare ghost-gluon vertices, the ghost triangle indeed 
yields a negative contribution leading to a zero crossing. The inclusion of the other diagrams, which 
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Figure 4. Comparison of different parameters for the three-gluon vertex model: Iiir = —1, A3 9 = 
0, 0.86, 1.08, 1.31, 1.53, 1.7 GeV. Top left: Symmetric configuration of the vertex model. Thick lines from left 
to right correspond to rising values of the parameter A3 9 . The thin line is the 1R leading contribution of the 
three-gluon vertex calculated from its DSE showing a zero crossing at very low momenta. Individual points are 
lattice data from Ref. [13]: Red circles correspond to iV = 16//3 = 2.2, green squares to N = 16//? = 2.5, blue 
diamonds to N — 22//3 = 2.2 and orange triangles to N — 22/ /3 = 2.5. Top right: Ghost dressing function. 
All curves lie almost on top of each other. Bottom: Gluon dressing function (left) and propagator (right). 
Curves from low to high correspond to values of A3 9 from low to high. 

are suppressed in the IR, may change its position slightly but they will not make it disappear. A 
more detailed investigation of the three-gluon vertex was done in Ref. [56] for two dimensions. There 
it was found that gluonic diagrams indeed do not affect the IR part. Note that the ghost propagator 
is hardly affected by the different three-gluon vertex models and consequently all curves for the ghost 
triangle contribution are essentially the same. 

5 Influence of a dynamical ghost-gluon vertex 

The ghost-gluon vertex is central for any truncation scheme of functional equations in the Landau 
gauge, where usually a bare vertex is employed. Deviations from tree-level might yield quantitative 
corrections in the mid-momentum regime, because here the difference to a bare vertex is largest. In 
this section we add the ghost-gluon vertex dynamically to the set of equations to be solved in order 
to assess how much quantitative change this entails for the propagators. For the three-gluon vertex 
we will use the model of Eq. (2.23). 
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Figure 5. Comparison of the results for the gluon propagator dressing function Z(p 2 ), the gluon propagator 
Z(p 2 )/p 2 and the ghost dressing function G(p 2 ) with lattice data [12]. The red/continuous line was obtained 
with a dynamic ghost-gluon vertex and the optimized effective three-gluon vertex, Eq. (2.23) with parameter 
set 2 of Tab. 1. For reference the green/dashed line is shown, which was obtained with the three-gluon vertex 
model of Eq. (2.17) and a bare ghost-gluon vertex. Lattice data is for f3 = 6 and lattice sizes of L = 32 (blue 
circles) and L — 48 (orange squares). 



5.1 Propagators 

In order to obtain a solution for the DSEs consisting of the propagator and the ghost-gluon vertex 
equations, the three-gluon vertex model plays a crucial role. For example, using the model of Eq. (2.17) 
it turns out that the gluon loop can dominate over the ghost loop in the gluon DSEs which results in 
a negative gluon dressing function. Thus the changes induced by the Bose symmetrization and the IR 
part are required in order to obtain a solution. Of course, the inclusion of two-loop diagrams might 
eventually change that. However, since they are neglected here, we use parameters for the three-gluon 
vertex model that yield the gluon propagator in agreement with lattice data as well as possible. We 
refer to this concrete vertex model in the following as optimized effective three-gluon vertex. 

We compare the resulting propagators to lattice results in Fig. 5. For reference we also show the 
result from a propagator only calculation with the three-gluon vertex of Eq. (2.17). Clearly, the gap 
between lattice and DSE results in the mid-momentum regime can be considerably diminished. The 
changes for the ghost dressing are small. For the gluon propagator the IR part looks overenhanced. 
However, this region is very sensitive to the scale setting, which was done here by placing the maxima 
of the gluon dressings at the same position as in the lattice data. 

Assessing the effect on the propagators of the ghost-gluon vertex alone is not straightforward, 
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parameter set 


hm 


A 3g 


ghost-gluon vertex 


1 


-1 


2.1 GeV" 


dynamic 


2 


-1 


1.8 GeV 


bare 


3 


-1 


2.1 GeV 


bare 



Table 1. Used parameters for the three-gluon vertex model. Sets 1 and 2 correspond to the optimized effective 
three-gluon vertex model for the respective choice of the ghost-gluon vertex. 

because the three-gluon vertex model depends on the propagators which are in turn changed by the 
ghost-gluon vertex. Thus the model is not the same even if the same parameters are used. The 
following comparison is therefore to be considered with a grain of salt. In Fig. 6 we show three 
solutions: The one with a dynamic ghost-gluon vertex and the optimized effective three-gluon vertex 
from above and two others with a bare ghost-gluon vertex. When using the same parameters with a 
bare ghost-gluon vertex, the mid-momentum enhancement of the gluon dressing function is too strong. 
Adjusting the parameters to the bare ghost-gluon vertex, however, shows that the three-gluon vertex 
model can incorporate even effects of the ghost-gluon vertex, as this solution then agrees very well 
with the one obtained from the dynamic ghost-gluon vertex. Furthermore one can also see that the 
ghost dressing function is almost unaffected by changing the three-gluon vertex parameters and that 
it is only slightly modified by the dynamic ghost-gluon vertex. The employed parameters for the 
three-gluon vertex model are summarized in Tab. 1. 

5.2 Ghost-gluon vertex 

Next we will present results for the ghost-gluon vertex. Up to now numerical non-perturbative in- 
vestigations of the ghost-gluon vertex were only done for specific momentum configurations [3, 13] 
or were restricted to the IR part [6]. In Ref. [9] a ghost-gluon vertex model based on an operator 
product expansion [8] was adapted such as to properly describe the ghost dressing via its DSE. Here 
we obtain the complete momentum dependence from a self-consistent calculation. Thereby it turns 
out that the structure is in general as simple as for special momentum configurations: The vertex is 
one for large momenta, develops a small peak in the intermediate momentum regime and then settles 
for a constant value in the IR. To illustrate this we show a few exemplary momentum configurations in 
Fig. 7. Additionally typical one-scale configurations are shown in Fig. 8. For the so-called asymmetric 
momentum configuration, where the gluon momentum vanishes, we compare our result with lattice 
data. More lattice results with better statistics exist for SU(2) [13]. However, in this SU{2) data the 
peak of the vertex dressing occurs at somewhat larger momenta than in our calculations. Whether this 
is a genuine difference between SU(2) and SU(3), whether it indicates an inconsistency in the scale 
setting or whether the effect might be due to our truncation employed for the vertex DSE remains to 
be further investigated. 

In addition to the ghost-gluon vertex obtained with the optimized effective three-gluon vertex 
model, for comparison also shown in Fig. 8 is a result obtained with the same hm but a considerably 
larger A3 9 . Both results yield a very similar ghost-gluon vertex. Therefore we conclude that the vertex 
is rather insensitive to the three-gluon vertex model in our truncation, although this could be expected 
as an indirect effect via the gluon propagator. If we proceeded one step further in the truncation of the 
ghost-gluon vertex DSE, by including the dressed three-gluon vertex, this will most likely change. If 
in a next step the three-gluon vertex is calculated dynamically this should yield a viable improvement 
of the ghost-gluon vertex truncation in the future. 
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Figure 6. Comparison of the results for the gluon propagator dressing function Z(p 2 ), the gluon propagator 
Z(p 2 )/p 2 and the ghost dressing function G(p 2 ). A dynamic (red/continuous line, parameter set 1) and a bare 
(green/dashed line, parameter set 2) ghost-gluon vertex are used together with the corresponding optimized 
effective three-gluon vertices. The blue/dotted line (parameter set 3) stems from a calculation using the same 
parameters as in the dynamic calculation and a bare ghost-gluon vertex. See Tab. 1 for details on the parameter 
sets. 

The system of equations can also be solved for scaling boundary conditions, viz. 1/(2(0) = 
[2]. By varying the boundary condition to a finite value a family of decoupling solutions can be 
obtained [5, 28, 50]. Above this value was chosen such to reproduce the lattice results for the ghost 
dressing function. Up to now it has not been investigated numerically if this ambiguity in the solutions 
also persists for the vertices. At least qualitatively a difference is expected [5]. However, the ghost- 
gluon vertex should become IR finite for both types of solutions. Any difference must therefore be 
of quantitative nature, fndeed we find such a difference as illustrated in Figs. 7 and 9 for selected 
momentum configurations: For scaling the IR value of the ghost-gluon vertex is larger than one, but 
for decoupling it is exactly one, i.e., all diagrams vanish in this limit. This can be understood from 
IR power counting of the dressing functions: In the Abelian diagram we have two ghost propagators 
and one gluon propagator. We get in the limit that all external momenta approach zero with one 
scale that the diagram scales as p 2 for decoupling, i.e., it vanishes, while for scaling the IR exponent is 
— 2k + 2k = 0. Thus the diagram is IR finite. The non- Abelian diagram with a bare three-gluon vertex 
is IR suppressed for both types of solutions. It is important to note that our result of a ghost-gluon 
vertex not approaching one in the IR is not in conflict with Taylor's theorem, since it only concerns 
the longitudinal dressing function as discussed in Section 2.2. 
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Figure 7. Selected momentum configurations of the ghost-gluon vertex for a decoupling solution. Fixed 
momentum (left) or angle (right) as indicated at the top of the plots. 

6 Summary and outlook 

In the present work we have extended typical truncation schemes of DSEs by including the ghost-gluon 
vertex dynamically together with the propagators in the calculations. The employed DSE truncations 
contain the IR and UV leading diagrams. The results for the vertex are in good agreement with 
available lattice data: Basically it stays flat with a small bump around 1 GeV. The influence on the 
propagators is most visible in the mid-momentum regime. As far as the different types of solutions for 
the propagators are concerned, this ambiguity is also present at the level of the vertex: For decoupling 
the vertex approaches one in the IR, while a scaling vertex receives IR contributions from the IR 
leading diagram thus settling at a value larger than one for asymptotically small momenta. 

The second improvement concerns the three-gluon vertex: While in many truncation schemes it 
is adjusted so as to obtain the correct UV anomalous dimension for the gluon propagator, we have 
devised a model that respects Bose symmetry and qualitatively captures the behavior seen in lattice 
calculations, especially the fact that its dressing function turns negative in the IR. Since only two- 
and three-dimensional lattice data is unequivocal in this respect, we have calculated the IR leading 
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Figure 8. The ghost-gluon vertex for one-scale momentum configurations. The red continuous line is the 
solution with the optimized effective three-gluon vertex (hm — —1 and A3 9 = 2.1 GeV), the green dashed 
one was obtained with the same hm but for A3 g = 2.6 GeV. Left: Symmetric momentum configuration 
p 2 = q 2 = k 2 . Right: The configuration with vanishing gluon momentum in comparison with lattice results 
from Ref. [12] (blue circles: N = 48 at @ = 6). 



contribution of the three-gluon vertex at the symmetric point. We did find a zero crossing, but at 
momenta lower than the lowest available momentum in lattice simulations. Thus the zero crossing 
itself is not relevant for the gluon DSE as all contributions in this regime get suppressed by the two 
attached gluon propagators. This is also corroborated by an analysis of the influence of the IR part in 
our model. However, the Bose symmetrization of the vertex leads to an increase of the mid-momentum 
strength of the gluon propagator. By choosing the parameters of the model appropriately, we could 
effectively incorporate the effects of two-loop diagrams in the mid-momentum regime. With this 
optimized effective three-gluon vertex we were able to close the gap to lattice data. 

An important next step in extending the truncation scheme is certainly the three-gluon vertex. If 
included self-consistently one can then disentangle contributions due to two-loop diagrams from the 
gluon loop. For the former the three-gluon vertex will also be required as it is contained in the so-called 
squint diagram. Furthermore, an extension of the truncation of the ghost-gluon vertex DSE could be 
to use the DSE with the bare vertex attached to the ghost leg. This then includes the three-gluon 
vertex as well. As an alternative to DSEs it might be interesting to look also at the ghost-gluon vertex 
equation of motion from an nPI effective action, where the bare vertices get dressed. This might 
compensate for the neglected diagrams in the present truncation. 

In summary, the presented extension of a standard truncation scheme confirms the qualitative 
reliability of previous schemes. Furthermore, it shows that by including higher vertex functions DSEs 
can successfully be used to obtain quantitative results, although the underlying system of equations 
is infinitely large. Thus DSEs provide a rather satisfying description of the non-perturbative regime 
which can be systematically improved. 
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Figure 9. Selected momentum configurations of the ghost-gluon vertex for the scaling solution. Fixed 
momentum (left) or angle (right) as indicated at the top of the plots. 
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A Kernels 

The kernel Kq of the ghost DSE given in Eq. (2.1) is 

Kc( P , q) = x2 + {y -f; 2x{y + z) (A.1) 
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with x = p 2 , y = q 2 and z = (p + q) 2 . The two kernels K 9 z h and K 9 J of the gluon DSE Eq. (2.2) read 
Kf{p,q) = - 



% 2 + (y — z ) 2 — 2s(y + z) 



K 9 J(p, q) = M J y2z2 (x 4 + 8-T 3 (y + z) + x 2 (-18y 2 - 32yz - I8z 2 ) + (A.2) 

(y - zf (y 2 + lOyz + z 2 ) + 2x(y + z) (Ay 2 - 20yz + 4z 2 ) ) . (A.3) 

The appearance of spurious quadratic divergences is discussed in Appendix B. The gist is that they 
are subtracted via additional terms that contain a damping function such the IR behavior is not 
influenced. For the respective damping scales in the ghost and gluon loops we used L\j V = 10 and 
2 in internal units for the ghost and gluon loops. The physical values can only be determined at the 
end of the calculation when the scale is set. 

Untruncated DSEs respect multiplicative renormalizability. The employed truncation and the 
method for subtracting the quadratic divergences do not change this as illustrated in Fig. 10. There 
we show the results for the dressings and the running coupling if different values for a(fj, 2 ) are chosen. 
This indirectly fixes the renormalization point /i 2 . For the comparison the values of the damping scale 
parameters Ljjv must have the same physical values, which can only be determined a posteriori. 

The kernels of the ghost-gluon vertex DSE Eq. (2.16) read: 

K Acc _ (k ■ rp ■ k — k 2 p ■ r)((p ■ r) 2 + k ■ r(—p 2 + p ■ r) — p 2 r 2 + p ■ k(p ■ r — r 2 )) 
1 ~ 2(p-r) 4 r 2 {k + r) 2 ' ( ' ' 



K? cc = 



i., — —((k ■ rp ■ k — k 2 p ■ r)((k ■ r) 2 p ■ r + r 2 (p ■ kp ■ r + (p ■ r) 2 — p 2 r 2 — p ■ kr 2 ) 

(p — r) 2 r 4 (k + rp \ 

+ k 2 (k ■ rp ■ r — (p 2 + p ■ k — p ■ r)r 2 ) + k ■ r(p ■ kp ■ r + (p ■ r) 2 — 2p 2 r 2 — 2p ■ kr 2 + p ■ r r 2 )) 

+ (fc 2 p 2 — p ■ k 2 )(k ■ r(2k ■ r p ■ r — k ■ r r 2 — p ■ k r 2 ) + k 2 (fc • r p ■ r + r 2 (—p ■ k — p ■ r + r 2 )))j . 

(A.5) 

They were generated automatically using the program DoFun [79, 80] and implemented in C++ with 
CrasyDSE [80]. The kernels given here and those for B(k;p,q) were compared with those given in 
Ref. [17]. Except for a minus sign in front of the Abelian diagram the expressions agree. We repeated 
the calculation of Ref. [17] with our kernels. The results are shown in Fig. 11. Note that the different 
sign leads to an IR value which also differs from that in Ref. [17]. 



B UV behavior of the propagator DSEs 

In the asymptotic UV regime the propagator dressing functions have a logarithmic behavior: 

G(aO=G(s)(wlog(~) (B.l) 
Z{x) = Z(s)(u\og{^) +l) 7 , (B.2) 

where 5 and 7 are the anomalous dimensions and s is some high UV momentum. If both one loop 
diagrams are taken into account, lu is HA r c a(s)/127r, 6 — —9/44 and 7 = —13/22 [60, 61]. We will 
show how these values emerge from the propagator DSEs with the model for the three-gluon vertex. 
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Figure 10. Test of multiplicative renormalizability, here for the scaling solution: The gluon and ghost dressing 
functions (top left and right, respectively) depend on the renormalization point, which is indirectly fixed by 
the choice of a(fj, 2 ). The red/continuous curve is for a(/j, 2 ) = 1, the green/dashed one for a{n 2 ) = 0.5. The 
running coupling (bottom) is renormalization group invariant and thus does not change. The small deviation 
in the UV is a numeric artifact: The physical values for Luv do not match exactly, because the scale is set at 
the end and we can only get approximately equal values. 




Figure 11. Dressing functions A(k;p, q) and B(k;p,q) after one iteration for the same momentum config- 
urations and the same analytic expressions for the propagators as in Refs. [3] from our kernels, i.e., with the 
signs of the kernels of the Abelian diagram changed. The main differences arise in the IR, where the Abelian 
diagram is dominant. 
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We follow the calculation outlined in Ref. [21], where the anomalous dimensions for the ghost-loop 
only truncation were calculated using the Brown-Pennington projector [81] for the gluon equation. For 
illustration purposes we use here as in Refs. [60, 61] a generalized projector: 

P^ip)=9^-(rf L - (B.3) 

At the end we set ( = 1 to recover the transverse projector. Acting with this projector onto the gluon 
DSE yields: 

— = Z 3 + N c g 2 J G(q 2 )G((p + qf)Kf^{p, q) 

+ n c9 2 [ z(q 2 )z(( P + q) 2 )K^(p,<i)r A3 (p,q,p + q) (b.4) 

Jq 

with the new kernels 



K ai,C ( \ = z% , ^(5x - x( + 4yQ x 2 (~19 + C) + 2xy(-17 + C)~18y 2 C 

z [P ' q> 24x 2 y 2 I2x 2 y 2 24x 2 y 2 

(x - y) 2 (x 2 + 10xy + y 2 () Ax 3 + xy 2 (-l7 + () + 4y 3 C - x 2 y(15 + C) 

2Ax 2 y 2 z 2 \2x 2 y 2 z 



(B.6) 



with x — p 2 , y — q 2 and z — (p + q 2 ). As also projects onto the longitudinal part of the ghost-gluon 
vertex, we already replaced it at this point by the bare expression. 

Here we are interested only in high external momenta p 2 , for which the integral is dominated by 
the region x < y < L 2 , where L is a UV cutoff. Since the dressing functions G{p 2 ) and Z(p 2 ) vary 
only slowly for high momenta, we replace (p + q) 2 by q 2 . This allows to calculate the angle integrals 
analytically: 

^r Z3+ '^f dyX -ir Lz(yW > (R7) 

1 v j. N <S* t\ -»K-4) + x«-2) „, ,„, , 
W) ~* Z ' + 195? J, * ~y G{y)G ^ 

384tH J x xy 2 

Here the Bose symmetric three-gluon vertex including the renormalization group improvement term 
from Eq. (2.23) but without IR part is used. The quadratically divergent parts are those proportional 
to ( — 4. We deal with them by adding terms to the kernels that cancel them. Since this is a UV 
problem and we do not want these terms to affect the IR or mid-momentum regime they are multiplied 
by appropriate damping factors. The kernels become 

K* h ><(p,q) -> K* h ><(p,q) ^( 4 " O^AnKv), (B.9) 

K^(p,q) K^(p,q) + i(4 - 0—fuv(y), (B.10) 

2 xy 
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where 

jW(y)=tanh(-^). (B.ll) 

Using these kernels every loop takes care of its own quadratic divergences. This approach is similar to 
the one adopted in Ref. [28], where the vertices were modified for this purpose. Since the quadratic 
divergences are only introduced due to the use of a gauge-variant regularization, they are completely 
artificial and can consequently be subtracted as done here. However, such a subtraction is not unique 
and can possibly change the finite parts. Thus we choose the parameter L in an interval for which the 
results depend least on it. 

Next we discuss the parameters a and (3. In order to reproduce the expected UV behavior, they 
have to fulfill the following condition: 

28 = 2a<5 + 7(2 + 2/3). (B.12) 

In order to fix them completely, we demand that this part of the vertex becomes finite in the IR, so 
that it does not interfere with the IR enhanced part, Eq. (2.21). i.e., 

(3 = (B.13) 

for decoupling and 

2(3 = a (B.14) 

for scaling. For the decoupling and scaling solutions we obtain then a = 3+1/8, (3 = and a = —2—68, 
(3 = — 1 — 38, respectively. 

The values for 8 and 7 can be obtained by performing the final integral after plugging in the 
perturbative expressions from eqs. (B.l) and (B.2) with s chosen such that G(s) = Z(s) = 1: 



(»o + r - * - (t)) ,+ ™ - (r I+7+ * 

(B.15) 

/ , /x\ \-r - ((C-2)-(24 + 0)A r c<7 2 17 , {L 2 \\ 1+2S ( , /.t\\1+2 
( wl0g W +1 ) = Z3+ 192^(1 + 28)l ( 1 + Wl0g b)j -( 1 + Wl ° g (l)) 



(B.16) 

The dependence on £ cancels between the two loops. The terms divergent for L — > 00 are canceled by 
the renormalization functions Z3 and Z3. From the rest the consistency conditions 

1+7 + 2(5 = 0, (B.17) 



64tt 2 (1 +7 + 8)oj 
13 N c g 2 
96tt 2 (1 + 2S)lu 

follow. We obtain for the anomalous dimensions and u: 



1, (B.18) 
1. (B.19) 



5—1, (B.20) 
7 = ~. (B.21) 



22 
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with a(s) = g 2 /4ir. 
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